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We study the computational complexity of the A'^-representability problem in quantum chemistry. 
We show that this problem is QMA-complete, which is the quantum generalization of NP-complete. 
Our proof uses a simple mapping from spin systems to fermionic systems, as well as a convex 
optimization technique that reduces the problem of finding ground states to A'^-representability. 
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The central theoretical problem in the field of many- 
body strongly correlated quantum systems is to find ef- 
ficient ways of simulating Schrodinger's equations: it is 
very easy to write down those equations, but notoriously 
difficult to solve them or even to find approximate solu- 
tions. The main difficulty is the fact that the dimension 
of the Hilbert space describing a system of N quantum 
particles scales exponentially in N. This makes a direct 
numerical simulation intractable: every time an extra 
particle is added to the system, the computational re- 
sources would have to be doubled. 

The situation is not hopeless, however, as in principle 
it could be that all physical wavefunctions, i.e., the ones 
that are realized in nature, have very special properties 
and can be parameterized in an efficient way. The idea 
would then be to propose a variational class of wavefunc- 
tions that capture the physics of the systems of interest, 
and then do an optimization over this restricted class. 
This approach has proven to be very successful, as wit- 
nessed by mean field theory and renormalization group 
methods. So far, an efficient variational class to describe 
complex wavefunctions such as those arising in quantum 
chemistry has not been found. 

One of the basic problems in quantum chemistry is to 
find the ground state of a Hamiltonian describing the 
many-body system of an atom or molecule. The essen- 
tial element that makes typical Hamiltonians very un- 
generic is the fact that at most 2-body interactions occur. 
This implies that the number of free parameters in such 
Hamiltonians scales at most quadratically in the number 
of particles or modes, and hence the ground states of all 
such systems form a small-dimensional manifold. 

In the case of a Hamiltonian with only 2-body interac- 
tions, the energy corresponding to a wavefunction is com- 
pletely determined by its 2-body correlation functions, 
and as a consequence the ground state will be the one 
with extremal 2-body reduced density operators. This 
fact was realized a long time ago, and led Coulson 1] to 
propose the following problem: given a set of N quantum 
particles, can we characterize the allowed sets of 2-body 



correlations or density operators between all pairs of par- 
ticles? 

If the particles under consideration are fermions, this 
has been called the N -representability problem (2|- Here, 
we consider the reduced density operators acting on pairs 
of fermions, and we want to decide whether they are con- 
sistent with some global state over N fermions. 

An efficient solution to the A'^-representability problem 
would be a huge breakthrough, as it would (for example) 
allow us to calculate the binding energies of all molecules. 
Therefore, a very large effort has been devoted to solving 
this problem and there exists a large literature on this 
subject (see e.g. HHIl). 

Here we will show that the iV-representability problem 
is intractable, as it is QMA-complete and hence NP-hard. 
The complexity class QMA (Quantum Merlin- Arthur) is 
the natural generalization of the class NP (nondetermin- 
istic polynomial time) to the setting of quantum comput- 
ing. Colloquially, a problem is in QMA if there exists an 
efficient quantum algorithm that, when given a possible 
solution to the problem, can verify whether it is correct; 
here the "solution" may be a quantum state on polyno- 
mially many qubits 23] . A problem is QMA- hard if it is 
at least as hard as any other problem in QMA; that is, 
given an efficient algorithm for this problem, one could 
solve every other problem in QMA efficiently fi^l- We 
say that a problem is QMA-complete if it is in QMA and 
it is also QMA-hard. 

In a seminal work, Kitaev JjJ proved that the Local 
Hamiltonian problem — determining the ground state en- 
ergy of a spin Hamiltonian that is a sum of 5-body terms 
(on n qubits), with accuracy ±e where e is inverse poly- 
nomial in n — is QMA-complete. In fact, it was later 
shown that this problem remains QMA-complete when 
restricted to 2-body interactions 0, 1^ , and even in the 
case of geometrically local interactions 

Because the Hamiltonians under consideration are lo- 
cal, the corresponding ground states have extremal local 
properties. The dual problem of determining the ground 
state energy of a local Hamiltonian is to decide whether 
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a given set of local density operators can be realized as 
the reduced density operators of the same global state. 
Checking the consistency of such a set of local density op- 
erators p^-^^\ i.e., checking whether there exists a global 
state a compatible with those local density operators, 
is a QMA-complete problem itself the present 

paper, we will use very similar techniques to prove that 
iV-representability, which is the fcrmionic version of that 
problem, is also QMA-complete. 

To sketch the main ideas of the proof, we will first con- 
sider the classical marginal problem: suppose we have N 
random variables that can take the values ±1, according 
to some joint probability distribution, and we define the 
2-variable marginal distributions 

p{ai,aj)^ ^ p(ai, 02, Cat). 

Given a set of marginal probability distributions, we 
would like to check whether they are consistent, i.e., 
whether there exists a global probability distribution 
whose marginals are equal to the ones we were given. 

Suppose that there exists an efhcient algorithm (whose 
running time is polynomial in N) to solve this problem. 
Then it would be possible to identify ground states of 
all Ising spin glasses. The strategy would be as follows: 
since the set of consistent marginals is convex and the 
energy is a linear function of the marginals, the problem 
amounts to minimizing a linear function subject to a set 
of convex constraints that can be checked in polynomial 
time. This can be done in polynomial time using the 
ellipsoid method ■ It therefore follows that an efficient 
algorithm for checking consistency allows to find ground 
states of Ising spin glasses, a problem that is known to be 
NP-hard 11^. Hence the problem of determining whether 
a set of binary marginals is consistent is itself NP-hard. 

Let us now formulate the iV-representability problem 
in the context of quantum chemistry. Electrons and nu- 
clei tend to arrange themselves such as to minimize their 
energy, and the binding energy of a molecule can be de- 
termined by calculating the minimal energy of the corre- 
sponding Hamiltonian. In practice, the nuclei are fairly 
well localized and can be treated as classical degrees of 
freedom, and the wavefunction of the N electrons can be 
approximated as a linear combination of tensor products 
of the d single-particle modes in the system (those form 
the basis set). 

As electrons are fermionic, the complete wavefunction 
must be antisymmetric, and this is most easily taken into 
account by working in the formalism of second quantiza- 
tion: 

1 

1^^)= E Q„...,,(alr...(aty^|17). 

n. ■■■.'<£ = 
n + ■■■ + M = « 

Here aj is the creation operator for the j'th mode, and 
|f2) represents the vacuum state without fermions. The 



creation and annihilation operators obey the following 
anticommutation relations: 

{ai, aj} = = {a|, a]} {oj, a]} = 2%. 

Note that we restrict ourselves to the subspace of states 
with exactly N fermions. d denotes the number of modes, 
which is typically much larger than N. The number of 
degrees of freedom is (^) , which grows exponentially in 
N when d > cN for some constant c > 1. 

In the case of quantum chemistry, the Hamiltonian 
typically contains only one- and two-body interactions 
between all modes, and so it can be written as a linear 
combination of terms of the form a,- and a]a\aiak- 

The 2-fermion reduced density matrix (2- RDM) is cal- 
culated by tracing out all but two of the fermions: 

p(')-tr3....,^pW. 

where p*-^-* is a mixture of states \il>) with exactly N 
fermions. In the language of second quantization, the 
matrix elements of the 2-RDM are given by: 

(2) I / t t \ 

The A'^-representability problem (with d modes) can 
now be stated as follows. Consider a system of N 
fermions and d modes, d < poly(iV). (For purposes of 
complexity, we consider N to be the "size" of the prob- 
lem.) We are given a 2-fermion density matrix p, of 
size '^^'^~^'> X ^^^^^-!^, where each entry is specified with 
poly(A^) bits of precision. In addition, we are given a 
real number /3 > 1/ po\y{N), specified with poly(iV) bits 
of precision. The problem is to distinguish between the 
following two cases: 

• There exists an A'^-fermion state a such that 
tr3,...,Ar(cr) = p- In this case, answer "YES." 

• For all A^-fermion states a, ||tr3....^jv(iT) — pHi > (3- 
In this case, answer "NO." 

If neither of these cases applies, then one may answer 
either "YES" or "NO." Note that we do not insist on 
solving the problem exactly; we allow an error tolerance 
of /? > l/poly(A^). We use the £i matrix norm or trace 
distance, ||A||i = tr to measure the distance between 
a and p. 

We will show that A^-representability is QMA- 
complete. The proof consists of two parts. First, we show 
that any 2-local Hamiltonian of spins can be simulated 
using a 2-local Hamiltonian of fermions with d = 2A^. 
Using techniques of convex programming, we show that 
an oracle for A^-representability would allow us to esti- 
mate the ground state energies of 2-local Hamiltonians; 
thus, A^-representability is QMA-hard. 

Second, we show that A^-representability is in QMA; 
specifically, we construct a quantum verifier that can 
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check whether a 2-particle state is iV-representable, given 
a suitable witness. 

Let us first show how to map a 2-local Hamiltonian, 
defined on a system of TV qubits, to a 2-local Hamiltonian 
on fcrmions, with d = 2N modes; this is the opposite of 
what has been done in 17} . The idea is to represent each 
qubit z as a single fermion that can be in two different 
modes ai,bi; so each A^-qubit basis state corresponds to 
the following A^-fermion state: 

i^i) ® • • • ® \zn) ^ ia\y-'^ {b\r ■ ■ ■ («jv)'^"" Ky" \^)- 

(1) 

Also, all the relevant singlc-qubit operators should corre- 
spond to bilinear functions of the creation and annihila- 
tion operators (this construction guarantees that opera- 
tors on different qubits commute). This can be achieved 
by the following mapping: 

< ^ alb, + bla,, <-> i (foja^ - a^b^) , < ^ 1 - 26|6,. 

(2) 

It can easily be checked that these fermionic operators 
obey the Pauli commutation relations within the sub- 
space spanned by states of the form (Q) . 

Recall that the qubit Hamiltonian can be written as a 
linear combination of Pauli operators. So the fermionic 
Hamiltonian can now be created by rewriting these Pauli 
operators with respect to the above mapping (note that 
only bilinear terms occur, and hence operators on differ- 
ent sites commute) . Notice that if the qubit Hamiltonian 
only contains 2-body terms, then the fermionic Hamilto- 
nian will only contain terms with at most 2 annihilation 
and 2 creation operators. 

The only thing that is left to do is to guarantee that 
there is exactly one fermion on every site i. This can 
be achieved by adding the following projectors as extra 
terms in the fermionic Hamiltonian: 

= (2a|a.- I)(26j6,- a). 

All the Pi are biquadratic and commute with all the op- 
erators introduced in Q , and hence the complete Hamil- 
tonian will be block diagonal. By making the weights of 
these projectors large enough (a constant times the norm 
of the Hamiltonian, which is at most polynomial in N), 
we can always guarantee that the ground state of the full 
Hamiltonian will have exactly one fermion per site |25j . 

Let us now assume that we have an efficient algorithm 
for iV-representability. We claim that this allows us to 
efficiently determine the ground state energy of a 2-local 
Hamiltonian on qubits, a problem which is known to be 
QMA-hard 0. We start by transforming the Hamilto- 
nian iJqubit on N qubits with 2-particle interactions into 
a Hamiltonian i/fcrmi on d = 2N fermionic modes with 
2-particle interactions (as described above). The prob- 
lem of finding the ground state energy of -fffcrmi can be 



expressed as a convex program, which can be solved in 
polynomial time using an oracle for 7V-representability 
(this is what we will show below). Our approach is similar 
to though in this case we are dealing with fermions 
rather than qubits. 

The basic idea is to construct a convex program 
that finds a 2-particle density matrix p which is N- 
representable, and which minimizes the expectation value 
of i?fcrmi- This program has polynomially many vari- 
ables, and it is easy to see that the set of 7V-representable 
states is convex, and (-fffcrmi) is a linear function of p. 

There are two minor complications. First, -fffcrmi de- 
scribes a system with 2N modes and an arbitrary number 
of particles, whereas A^-representability pertains to a sys- 
tem with exactly N particles. But this is not a problem, 
because the interesting behavior in iJtcrmi occurs in a sub- 
space where all the states have exactly N particles. So, 
in our convex program, we only need to consider states 
p that have exactly N particles. 

This lets us simplify -fffcrmi as follows. Since we are 
only interested in how it acts on A'^-particle states, we 
can write ajaj = jfrrjaliJ^k we can assume 

that all the terms in -fffcrmi are of the form a|a]a;afe. 

Second, convex optimization algorithms usually re- 
quire that the set K of feasible solutions be full- 
dimensional, i.e., K cannot lie in a lower-dimensional 
subspace. So we have to represent the state p in such 
a way that there are no redundant variables. 

Let 5 be a complete set of 2-particle observables jl^ , 
and let £ = \S\; note that e < poly(d) < poly(iV). We 
represent p in terms of its expectation values as = tr(S'p) 
for all S E S; let a E denote the vector of these 
expectation values. We define K to be the set of all a 
such that the corresponding state p is iV-representable. 
Note that the 7V-representability oracle allows us to test 
whether a given point cf is in if . 

We write our Hamiltonian in the form H — X)sg5 "^^S 
(plus a constant term); the coefficients 75 can be com- 
puted using the Gram-Schmidt procedure. It is easy 
to see that tv{Hp) — TliSeS^sois- Our convex pro- 
gram is as follows: find some a E K that minimizes 
/(a) = Eses^s^s- 

We can solve this convex program in polynomial time, 
using the shallow-cut ellipsoid algorithm (see 0], and 
references cited therein). We mention a few technical 
details. The algorithm requires some additional informa- 
tion about K , namely a guarantee that K is contained 
in a ball of radius R centered at 0, and K contains a ball 
of radius r centered at some point p. (Also, the running 
time of the algorithm grows polynomially in log(i?/r).) 
In our case, we can set R = and r — 1/ poly(i?) [27|. 

There is also the issue of numerical precision. Our or- 
acle for A^-representability has limited accuracy — it may 
give incorrect answers for points near the boundary of 
K . On the other hand, we only need an approximate 
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solution to the convex program, in order to solve the Lo- 
cal Hamiltonian problem. It turns out that 1/ po\y{N) 
precision is sufficient. The ellipsoid algorithm still works 
in this setting; see for details. 

Note that, in place of the ellipsoid algorithm, we could 
have used a different algorithm based on random walks 
in convex bodies [3 i this was the approach used in |0| . 
However, it is not clear if we could use one of the interior- 
point methods for convex optimization; these methods 
are substantially faster, but they usually require an ex- 
plicit description of the constraints, not just a member- 
ship oracle. 

This completes the proof that A^-representability is 
QMA-hard. As a corollary, we have also proven that 
estimating the ground state energy for local fermionic 
Hamiltonians is QMA-hard. Note that the use of the el- 
lipsoid algorithm to reduce a convex optimization prob- 
lem to a convex membership problem is not new; see 
for references to related work. 

Next, we show that iV-representability is in QMA. 
That is, we construct a poly-time quantum verifier V 
that takes two inputs: a description of the problem (that 
is, p and /?); and a "witness" t, which is a quantum state 
on polynomially many qubits. The verifier V should have 
the following property: if p is iV-representablc, there ex- 
ists a witness t that causes V to output "true" with 
probability > pi', if p is not A^-representable (within er- 
ror tolerance f3), then for all possible states t, V outputs 
"true" with probability < po; and Pi — Po > 1/ poly(-/V). 

The idea is that, when p is iV-representable, the correct 
witness r consists of (multiple copies of) an iV-fermion 
state a that satisfies tr3^...^jv(o') = p. Then the verifier 
can use quantum state tomography to compare a and p. 

We represent the iV-fermion state a using d qubits, via 
the following mapping: 

Call the resulting qubit state a. We use the Jordan- 
Wigner transform to map the fermionic annihilation op- 
erators to qubit operators: 

a, ^ A, = -{®k<t CTfe)® |0)(lk- 

Thus, an observable O = ajajajafe -t- a^ajajUi is trans- 
formed into 6 = A\A^^AiAk AlAjAjAi, which is a 
tensor product of many single-qubit observables and one 
four-qubit observable. 

We claim that the expectation value (O) can be es- 
timated efficiently. Without loss of generality, assume 
that the eigenvalues of O lie in the interval [0, 1]. Then 
there is an efficiently-implementable measurement which 
outputs "1" with probability (O) j2SiJ. By repeating this 
measurement on multiple copies of the same state, we 
can estimate (O); for our purposes, it is enough to have 
polynomially many copies of the state. 



We now describe the verifier V. The witness r consists 
of several (i.e., polynomially many) blocks, where each 
block has d qubits, supposedly representing one copy of 
the state a. On each block, V measures the observable 
J2k and if the outcome does not equal N, V out- 

puts "false." This projects each block onto the space of 
iV-particle states. 

Next, V performs measurements on each block, to es- 
timate the expectation values of a, for a suitable set of 
observables. Then V checks whether they match the ex- 
pectation values of p. One problem arises: the prover 
could try to cheat by entangling the different blocks of 
qubits. One can show that this does not fool the veri- 
fier, using a Markov argument, as was done in This 
suffices to show that iV-representability is in QMA, and 
hence finishes the proof. 

What can be said about the complexity of the pure- 
state A'^-representability-problem, where one has to de- 
cide whether the reduced density operators arise from 
a pure state of N fermions? In that case, the ver- 
ifier must be able to convince himself that the state 
he gets is pure. This can be done when he gets two 
states p and a that are promised to be uncorrelated, 
i.e. that he gets the state p ® a: then Arthur can cal- 
culate the expectation value of the observable tr(pCT), 
and this can only be close to 1 when p is pure and 
a c:^ p. Indeed, if tr((T^) < 1 — £, then for all states 
T, tr(crT) < ^tr(cr2)tr(T2) < (1 - ef/^ < 1 - e/2. The 
problem however is that Merlin can cheat, and hand out 
a correlated state to Arthur, such that the above test 
becomes inconclusive. This is precisely the feature that 
distinguishes the complexity class QMA(k) from QMA 
fie]: in QMA(k), the verifier is promised to get a tensor 
product of different states, and it has been conjectured 
that QMA(k) is strictly larger than QMA. The above 
discussion shows that the pure-state A^-representability 
problem is contained in QMA(k), and it is hence plausi- 
ble that it is harder than the A^-representability problem. 
It would be very interesting to investigate whether the 
problem is QMA(k)-complete. 

It is remarkable that checking consistency of 2-body 
reduced density operators of fermionic states is so hard, 
while checking consistency of 1-body reduced density op- 
erators is simple 2] . This can be easily understood from 
the previous discussion: the extreme points of the con- 
vex set of 1-body density operators {a\aj) correspond 
to ground states of Hamiltonians only containing bilin- 
ear terms in a\ and a^; such Hamiltonians can easily be 
diagonalized as they represent systems of free fermions, 
and hence the consistency problem can easily be solved. 
As shown in 0, consistency can be decided in that case 
based solely on the eigenvalues of the reduced density 
operators. A number of related problems, where con- 
sistency only depends on the eigenvalues, have been in- 
vestigated recently Il7l| ; most remarkably, a characteri- 



5 



sation has been obtained for the polytope of 1-particle 
marginals that are iV-representable under the condition 
that the iV-particle wavefunction must be pure . 

These results have to be contrasted with our problem 
of finding the A^-representable 2-body density operators, 
where the eigenvalues alone are not enough to decide 
consistency but also the eigenvectors are relevant. Ac- 
tually, let us consider the simpler problem of deciding 
A^-representability of 2-fermion density operators where 
only the diagonal elements Dij = {a\a}-ajai) are speci- 
fied. If we consider the case d = 2N and the mapping 
discussed above, one easily finds that the extreme points 
of this set would be obtained by ground states of local 
spin Hamiltonians which only contain commuting op- 
erators. These correspond to spin-glasses, and so the 
problem of deciding A^-representability of {Dij} is NP- 
hard [l^ . It was indeed pointed out a long time ago that 
7V-representability restricted to the diagonal elements is 
equivalent to a combinatorial problem 19] that was later 
shown to be equivalent to the NP-hard problem of decid- 
ing membership in the boolean quadric polytope H^l • 

Let's finally discuss the relevance of the above results in 
the context of quantum chemistry. We have shown that 
it is a hopeless task to determine the ground states of 
all local fermionic Hamiltonians, and in particular, that 
an approach by means of the A'^-representability prob- 
lem is intractable, even on a quantum computer. But 
it is possible that other physical systems, e.g. systems 
with different particle statistics, additional symmetry or 
a limited number of modes, might allow for an efficient 
characterisation of the two-particle reduced density ma- 
trices, and hence an efficient calculation of the ground 
states of local Hamiltonians. This seems to be the case 
in e.g. one-dimensional translational invariant spin sys- 
tems, where the density matrix renormalization group 
[2l| allows for a systematic approximation of the allowed 
convex set of reduced density operators from within |22l | . 
A very different numerical approach has been developed 
in the context of quantum chemistry and it known as the 
the contracted Schrodinger equation Q. In that method, 
one approximates the convex set from the outside, and 
the 7V-representability problem is a crucial ingredient of 
the algorithm. The foregoing discussion shows that this 
approach has to break down in the most general case, and 
it would be very interesting to investigate the conditions 
under which those approximations are justified. 

In conclusion, we investigated the problem of A^- 
representability, and characterized its computational 
complexity by showing that it is QMA-complete. Ob- 
viously, the theory of quantum computing was a prereq- 
uisite to pinpoint the computational complexity of this 
classic problem. . 
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following observables: Xij = a\aj + a}jai, for all I -< J; 
Yij = —ia\aj + ia'jai, for all I < J; and Zi — a\ai, for 
all I except the last one. These operators are Hermitian, 
with eigenvalues in the interval [—1, 1]. Taking real linear 
combinations, these operators form a basis for the space 
of 2-fcrmion density matrices. 
[27] It is easy to see that K is contained in a ball of radius 
R = since for all a £ K, we have — 1 < as < 1. 

The second claim, that K contains a ball of radius 
r = 1/ poly(/), is loss trivial. The proof is as follows. 

We will consider 7V-representability for different values 
of N; let Kn denote the set of all vectors d that are 
A'^-representable. Obviously, K2 contains a ball of radius 
1/ poly(^) (this is the trivial case). We will prove that Kn 
contains a ball of radius 1/ poly(£), for all 2 < A'^ < d — 2. 

We define "particle-hole" observables, by replacing ai 
with a|, and vice versa: X'jj = a/Oj + aja]; Y^j = 
—iaia\ + iaja\; and Z'j = aia\. Let a denote a vector 
containing expectation values for these observables; let 
K'pj be the sot of all a' that are A^-rcpresentable. 

First, we claim that K2 ~ K'^_2- Notice that there is a 
natural correspondence between the 2-particle Slater ba- 
sis states and the (d — 2)-particle Slater basis states: the 
2-particle state with modes i and j occupied corresponds 
to the (d — 2)-particlc state with modes i and j empty. 

So take any point a £ K2, which represents the ex- 
pectation values of the 2-particle observables for some 2- 
paxticle state a. Use a to construct a (d— 2)-particle state 
r, by replacing each 2-particle Slater basis state with the 
corresponding (d— 2)-particle Slater basis state. Then the 
expectation values of the 2-particle observables for cr are 
exactly the expectation values of the 2-hole observables 
for T. So a is in This shows that K2 C K'^_2- A 

similar argument shows that Q K2. So K2 = K'j^-2- 

Next, we claim that there is an invertible linear trans- 
formation A that maps K'^_2 to Kd-2- 

First we map K'^_2 to Kd-2- Observe that we can 
write each 2-particle operator as a linear combination 
of 2-hole operators. (This holds provided we restrict 



the operators to act only on the subspace of (d — 2)- 
particlc states.) For instance, when / n J = 0, a|aj = 
aja\. When |7 n Jj = 1, we write equations such as 
alajaiaj = alajajai — — ajaj(l — aiol) = —aja\ + 
aja\aia\ = —aja\ + aja^ajaj; then use the fact that 
o-jil = {llk^i,ja'kal)aja\ = Efc^^ij afc%ai4- When 
/ = J, we have a\ai — X]i:,n/=0 o.lcl\. (For each of these 
equations, it is easy to check that the left and right sides 
act identically on all (d — 2)-particle Slater basis states.) 

Thus, for any (d — 2)-particle state a, the expectation 
values of the 2-particle observables are linear functions of 
the expectation values of the 2-holc observables. Thus we 
have a linear transformation A that maps K'^_2 to Kd-2- 

Similarly, we can map Kd-2 to K'd-2- Wc write each 
2-hole operator as a linear combination of 2-particle oper- 
ators (again, restricting the operators to act only on the 
subspace of (d— 2)-particle states). For instance, when in 
J = 0, aio^j — a]jai. When |/n Jj = 1, we write equations 
such as ajaiojaj = aiO^^aia] = — ojai(l — o|o;) = — ojoj-l- 
Ojaiajai — —a'jai + ajaj^ajOi; then use the fact that 

= i-d^Ek^i,jalak)a]ai = ^ Efc^,,, 
When I = J, we write aia\ = ai-ia\_^ai2a\^ = (1 — 
"ii«n)(l ~ ot^tti2) = 1 - ot^«ii - al2«i2 + OjOj; then 
use the fact that 1 — ^2^) El o-^il-, and a\ai = 
dTs Ei^i a\o'i4<iv = -dh T.ijLi ajajaitti. 

Thus, for any (d — 2)-particle state a, the expectation 
values of the 2-hole observables are linear functions of the 
expectation values of the 2-particle observables. Thus we 
have a linear transformation B that maps Kd-2 to K'd_2, 
and B = ^"^ 

We claim that Kd-2 contains a ball of radius 1/ po\y{£). 
We know that K'^-2 contains a ball of radius l/poly(£) 
(since K2 = ^^^-2)1 and we will show that the map A 
docs not shrink this too much. Write the singular value 
decomposition A = UDV, where U and V are unitary, 
and D is diagonal, with diagonal entries Da > 0. Let B = 
A~^. Looking at the matrix elements of B, we can see 
that tr(BtS) = Eij \Bij\'^ < poly(^). At the same time, 
tr(B+B) = tr{UD-WV-'^D-'^U-'^) = tr{D-^) > Dj:^, 
for all i. So we have Da > 1/ poly(i!), for all i. Thus the 
map A shrinks the radius of the ball by at most a poly(£) 
factor. So Kd-2 contains a ball of radius l/poly(£). 

Now we are almost done. Notice that Kn-i 5 Kn, 
for all 3 < A'^ < d — 2 (since, if a is consistent with 
some A^-fcrmion state cr, then it is consistent with the 
(A — l)-fcrmion state a' that results from tracing out 
the last particle). So we conclude that Kn contains a 
ball of radius l/poly(£), for all 2 < AT < d - 2. This 
completes the proof. 
[28] For example, let Ai and \9i} be the eigenvalues and eigen- 
vectors of O, and note that they arc easy to compute. 
Add an ancilla qubit, perform the unitary operation 
U : \ei)\0) i-> ]6»i) (VI- Ai]0) -F\/A7|1», and measure 
the ancilla in the 0,1 basis. 



